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Abstract 

Multiscaling properties of the Anderson localization of the cosmic 
electromagnetic fields before the recombination time are studied and 
results of a numerical simulation for a random banded matrix ensemble 
are found to be in good agreement with the large-scale cosmic data. 
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1 Introduction 



Large-scale electromagnetic fields are expected to be generated by the cos- 
mological phase transitions (electroweak and QCD) [T|-[T2| before the recom- 
bination time. On the other hand, before the recombination time the cosmic 
plasma is strongly fluctuating at the large scales due to the gravitational 
density perturbation. These plasma density fluctuations can play a role of 
randomly distributed scatterers for the electromagnetic fields and make it 
impossible for the electromagnetic fields to propagate if the wave length is 
longer than the elastic mean free path of electromagnetic wave by the scatter- 
ers. The localized modes are thermally excited and the electric field is more 
or less aligned in the localization region. The localization phenomenon by 
random fields was first proposed by Anderson to explain the behavior of 
the electric resistivity of metals in the presence of impurities and the theories 
were developed by himself and his collaborators [H],[ini- The localization 
phenomena were then observed in various physical systems. Probably the 
first paper which studied the possible localization of electromagnetic field in 
plasma is the Ref. JIHJ (see also [T7|). Roughly speaking, the Anderson lo- 
calization phenomenon occurs because of the destructive interference of the 
incident wave and the scattered waves from the randomly distributed sources. 
It is known that the localization occurs if the strength of the randomness is 
larger than a critical value. For a sufficiently large random potential, only 
the modes of the wave length shorter than the mean free path can propa- 
gate, while the other modes are localized. The authors of Ref. jTH] argue 
that this localization phenomenon should occur in the matter dominated era 
(before the recombination time) in the standard cosmology. The authors of 
[TH] took into account a single (Jeans) scale l c for the primordial perturbation 
for simplicity. Inclusion of a scale invariant spectrum for scales larger than 
l c produces a hierarchy of structures, which needs in a multifractal descrip- 
tion. During the past decade, multifractality of critical eigenfunctions at the 
Anderson transition has been a subject of intensive analytical and numerical 
studies (see Refs. jTHJ-EZl and references therein). A remarkable virtue of the 
Anderson localization mechanism is that the photon diffusion does not occur 
to the localized mode. This makes the localized structures to be imprinted 
in the photon gas after the decoupling of the photon gas (so-called cosmic 
microwave background, or CMB) from the baryonic matter. In the present 
paper we will show that results of the numerical simulations of multifractal 
Anderson localization are consistent with the CMB data. 
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2 Anderson localization of cosmic electromag- 
netic fields 



Let us, following to the paper ^H]; consider Anderson localization of cos- 
mic electromagnetic fields for the spatially flat Friedmann-Robertson- Walker 
Universe. We assume that the charged particles are non-relativistic. So the 
four velocity of the charged particles: u a = (a" 1 ,!;*), with v l being small 
compared with unity. Equation for the vector potential A in terms of the 
conformal comoving coordinates (A,x) in this case is 

d 2 A d 2 A _ e 2 na 2 A 
d\ 2 <9x 2 m 

where A is the conformal time defined by A = / — , m is the electron mass 
and e is the charge, n is the density. 

Split the plasma density n into the homogeneous part n and the space 
dependent random part 5n; 

n = n + Sn (2) 
Then the right hand side of the plasma equation becomes 



[^ + ^^]A (3) 
a an 



where u) r ^ c = y e ^ rec . Here we have used the conservation law: n = n rec /a 3 
before the recombination time, with the scale factor a normalized at the re- 



combination time. n rec is the electron density (therefore the baryon density) 
at the recombination time. The first term in the bracket is spatially constant 
but a decreasing function of the conformal time A. The second term repre- 
sents the fluctuating random environment induced by the baryonic density 
perturbation ^ which is equal to = in our non-relativistic case. Knowing the 

standard theory of linear perturbation in cosmology we see that = = ^ oc a 
in the matter dominated era, so that the second term is independent of the 

conformal time in that era. We write the random potential as V(x) = p a -=r, 
which is a function of comoving coordinates x only. So far we have discussed 



3 



the period after the equal time 1 and before the recombination time and found 
that the random potential induced by the density perturbation is constant 
in time. As we shall see this is important for the stability of the Anderson 
localization. For the matter dominated universe, the scale factor is given by 
a(t) = (tM 2/3 = (^M 2 with A = / -%r = t?£3fV3. One can reduce the 

v ' v tree ' y Strec' J a(t) rec 

Maxwell equations to the following wave equation as 

<9 2 A 8 2 A 

+ H = m * + „(*)]A (4) 

Here Q 2 = 9tl ec u>p ec2 . We expand the vector potential A in terms of the 
complete ortho normal set {ip n } as 

A = S a n/n(A)^„(x) + C.C. (5) 

with a„ being arbitrary coefficients. Here ^„(x) is a normalized solution of 
the eigenvalue equation: 

[-A + V(x)]^x) = £> n (x) (6) 

Localization means that a number of modes {ip n ,n < ^o} are bound 
states exponentially falling off at infinity. The marginal state no corresponds 
to the so-called mobility edge as we shall explain further. 



Here we only consider the adiabatic baryonic density perturbation for 
simplicity, the spectrum of which is roughly like the picture below. To sim- 
plify our discussion we assume in this section that the spectrum is sharply 
peaked around the Jeans length l c at the equal time and the distribution of 
perturbation is Gaussian. Inclusion of the scale invariant spectrum for scales 
larger than l c produces a hierarchy of structures (see next section). 

To understand the Anderson localization in the Universe before the re- 
combination time, we have to evaluate the mobility edge of the wave number, 
which is roughly the inverse of the mean free path l m f of the electromagnetic 
wave by the random potential. Consider a Schrodinger equation for the spa- 
tial part of the electromagnetic field: 

[-A + V(x)]ij> = E 2 tp (7) 

1 The equal time t eq is the time when the matter and radiation energy density are 
equal. This value depends on the value of the present baryon density. We take = 
(h^c )2/3 _ g ag a typical value for Q — 0.2. 
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with x being the comoving coordinates. Recall that the random potential 
V(x) = ujp ec2 = does not depend on the conformal time in the matter dom- 
inated era. The characteristic length scale l c of the density perturbation is 
given by the Jeans length at the equal time, which is roughly the same as 
the horizon scale at that time. After that time the length scale l c develops 
as oc a. This gives l c ~ 10 23 cm at the recombination time, while the horizon 
is ~ 2 x 10 23 cm. Then the mean free path is given by l m f = — where a 
is the elastic scattering cross section of the electromagnetic wave by a single 
potential and n c ~ ^r^ 3 is the density of "impurities". The elastic cross sec- 
tion is estimated roughly as ~ 4ti1 2 by the standard wave mechanics, because 
the "energy" is low and is roughly given by E 2 l 2 pa 1, while the potential 
is effectively gigantic V x I 2 pa uj^ 2 ^l 2 pa 5 x 10 32 . Here we have taken 
— ~ 10~ 5 because the assumed baryonic perturbation can be comparable to 
the temperature fluctuation of CMB observed by COBE (see [2E| and next 
section). However, the localization length and other results in the present 
consideration are insensitive to the magnitude of the density perturbation. 
Therefore the localization should occur, since the critical value of VI 2 will 
be of order one. Then the mean free path is given by l m f ~ l c /3 and the 
mobility edge is E* ~ which is consistent with our low energy picture. 

Fortunately in the present case, the randomness is huge so that an esti- 
mate of the localization length is available |30]. We have £ = l c j log( | / i? 2 ) , 
where l c is the minimum length scale in which range the potential can be 
considered as a constant. For the adiabatic baryonic perturbation, it will be 
reasonable to take the Jeans length at the equal time t eq as the small scale 
cut-off l c , which is ~ 10 23 cm at the recombination time. In the present case, 
|l/|/i? 2 pa 5 x 10 32 so that £ ~ 10 21 cm at the recombination time. 

One can estimate the magnitude of the thermally excited localized elec- 
tric fields as E 2 oc pa T/£ 3 in a localization region of the size £. Here the 
contribution from a single localized mode only has been taken into account, 
because the peaks of the other localized modes are randomly distributed and 
will be located outside of the region under consideration. The estimate of 
E 2 oc can be derived, for instance, by the equal partition law. Note that the 
propagating modes(n > no) contribute to the Stefan-Boltzmann law. At this 
stage we realize that approximately T/u p photons occupies a localized state. 
As far as T/uj p » 1, a coherent localized state of the size £ is a good picture 
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so that we can approximately describe it by a classical field: 

no 

Eio C = " E V T M/n(^H(x)^n(x) + C.C. (8) 

where e n (x) is the polarization vector, which is perpendicular to the gradient 
of the wave function ip n (x). Therefore, in a localization domain of the size 
£ in the cosmic plasma, we have aligned electric and magnetic fields, which 
are perpendicular to each other and oscillate almost at the plasma frequency. 
One might worry about the damping of the electromagnetic structure of the 
size £ by photon diffusion. However, the diffusion does not occur to the 
localized modes. This is one of the virtues of the Anderson localization 
mechanism. Since the localization is a delicate interference phenomenon, the 
random potential V(x) should vary slowly so that the motion of wave packet 
is not disturbed ^H]- This condition is met before the recombination time. 
In the previous consideration we, following to [TH], only took into account a 
single scale l c for the primordial perturbation for simplicity. The contribution 
from the larger scales to the localization phenomenon is interesting in the 
sense that it may account for the hierarchy of the structures. 

3 Multifractal Anderson localization 

Even arbitrary close to the mobility edge a state should occupy only an 
infinitesimal fraction of space if it is to be labeled a localized state. On the 
other side of the mobility edge the states should extend through the whole 
volume. Both characteristics can be accommodated at the mobility edge if 
one assumes a fractal wave function with filamentary structure likes a net 
over the whole volume, as suggested originally by Aoki [2H]. Observations of 
the scaling properties of different characteristics of the eigenstates show that 
the wave functions cannot be adequately treated as a simple fractal. Rather, 
the more general concept of multifractality has to be employed, yielding a 
set of generalized dimensions. 

The multifractal fluctuations of eigenf unctions can be characterized by a 
set of inverse participation rations 

P p = J d d r\^(r)\ 2 ? (9) 

with an anomalous scaling with respect to system size L 

P p ~L^, r p = D p (p-l) (10) 
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The scaling (10) characterized by an infinite set of generalized dimensions D p 
implies that the critical eigenfunction represents a multifractal distribution. 
During the past decade, multifractality of critical eigenfunctions at the An- 
derson transition has been a subject of intensive numerical studies see Refs. 
JH|-[27] and references therein. In recent paper |2Zj the authors explore the 
fluctuations at criticality for the power-law random banded matrix (PRBM) 
ensemble. The model is defined as the ensemble of random Hermitian NxN 
matrices. The matrix elements are independently distributed Gaussian 
variables with zero mean (Hy) = and the variance 



W) = A\i-j)\), (11) 



where a(r) is given by 



a 2 (r) = , — - (12) 

At a — 1 the model undergoes an Anderson transition from the localized 
(a > 1) to the delocalized (a < 1) phase. At a = 1 the PRBM model 
was found to be critical for arbitrary value of b; it shows all the key fea- 
tures of the Anderson critical point, including multifractality of eigenfunc- 
tions. In a straightforward interpretation, the PRBM model describes an 
one-dimensional sample with random long-range hopping, the hopping am- 
plitude decaying as l/r a with the length of the hop. Such a random matrix 
ensemble arises in various contexts in the theory of quantum chaos |33j . [3*1] 
and disordered systems [SSI , EH] , jSZI • l n the P a P er the authors studied 
the periodic (critical) generalization of (12) 

a -2(r) - i | 1 *™ 2 (WAQ fl3l 
{T) - L+ V (tt/AO 2 [ } 

Indeed, for iV ^> r Eq. (13) reduces to Eq. (12) at the critical case. The 
generalization (13) allows diminish finite-size effects (an analog of periodic 
boundary conditions). At b = 1 a limiting distribution can be reached most 
easily and, therefore, the finite-size effects are minimal |27j . This value of 
parameter b corresponds to crossover from weak (b ^> 1) to strong (1 > 6) 
multifractality. Below we will compare results of the numerical simulations 
performed for b = 1 with the large-scale CMB data. 
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Temperature (T) dissipation rate of the photon (CMB) gas can be char- 
acterized by a "gradient" space measure (jHE], P- 381): 

* - ^2!* (14 ) 

where v r is a subvolume with space-scale r. Scaling law of this measure 
moments, 

(y p ) 

^4- ~ ^ (15) 

(Xr) P 

(where (...) means an average) is an important characteristic of the dissi- 
pation rate [HEMHS]- The exponents be related to correspondent 
generalized dimensions by equation [39J: 

/ip = (d-D p )(p-l) (16) 

The COBE satellite was launched on 1989 into 900 km altitude sun- 
synchronous orbit. The Differential Microwave Radiometers (DMR) oper- 
ated for four years of the COBE mission and mapped the full sky. The 
instrument consists of six differential microwave radiometers, two nearly in- 
dependent channels that operate at each of three frequencies: 31.5, 53 and 
90 GHz. Each differential radiometer measures the difference in power re- 
ceived from two directions in the sky separated by 60°, using a pair of horn 
antennas. Each antenna has a 7° beam. In this investigation we use a DMR 
data map of the cosmic microwave radiation temperature. The COBE-DMR 
team made an effort to remove the most of the Galactic microwave emission 
(i.e. thermal emission of cosmic dust, free-free emission and synchrotron 
emission) from the map using so-called "combination technique" in which a 
linear combination of the DMR maps is used to cancel the Galactic emis- 
sion. This was made in order to extract the cosmic microwave background 
radiation. 

To organize the DMR data in the temperature maps the sky was divided 
into 6144 equal area pixels. These are formed by constructing a cube with 
each face divided into 32 x 32 = 1024 squares, projected onto a celestial 
sphere in elliptic coordinates. The projection is adjusted to form equal area 
pixels having a solid angle of 47r/6144 sr or 6.7 square degrees. Because the 
7° beamwidth of the sky horn is greater than separation between pixels (2.6° 
average), this binning oversamples the sky. 
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Using the CMB pixel map, we will calculate the CMB temperature gra- 
dient measure using summation over pixel sets instead of integration over 
subvolumes v r . So that multiscaling (15) (if exists) will be written as 

Us/ ^ s -mp (17) 



(Xs) 



where metric scale r (from (15)) is replaced by number of the pixels, s, char- 
acterizing size of the summation set. The Xs is a one-dimensional surrogate 
of the real dissipation rate \ r . It is believed that the surrogates can repro- 
duce quantitative multiscaling properties of the dissipation rate j^Hl- Since in 
our case (x s ) is independent on s, we will calculate the exponents fi p directly 
from scaling of (x?)- 

Figure 1 shows scaling of the CMB temperature dissipation rate moments 
(x p s ) calculated for the DMR map. The straight lines (the best fit) are drawn 
to indicate the scaling in the log-log scales. 

Figure 2 shows the generalized dimensions D p extracted from figure 1 
(circles) using Eq. (16). Crosses in this figure correspond to the general- 
ized dimensions calculated for the above described PRBM model j2Zj of the 
Anderson localization. 



4 Discussion 

The Anderson localization is an unique linear phenomenon which exhibits 
profound multifractality. This allows explain the observed multifractality 
of the CMB dissipation in the frames of standard (linear) approach to the 
CMB, which seems to be the most compatible with the small magnitude of 
the CMB fluctuations. Generation of the cosmic electromagnetic fields in the 
electroweak or QCD phase transitions occurring in the particle physics scales 
needs in a non-linear (turbulent) mechanism for reasonable space scale expan- 
sion of the fields pQ- [T2]. Then the Anderson localization of the large-scale 
electromagnetic fields can be considered as an effective linear mechanism of 
the final CMB fractalization. This is because the photon diffusion, which can 
effectively suppress the non-linear fluctuations just before the recombination 
time, does not occur to the localized modes. 
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Figure Captions 

Figure 1. Scaling of the CMB temperature dissipation rate moments (x^) 
calculated for the four-year DMR map. The straight lines (the best fit) are 
drawn to indicate the scaling in the log-log scales. 

Figure 2. Generalized dimensions D p extracted from figure 1 (circles) 
using Eq. (16). The crosses correspond to the generalized dimensions calcu- 
lated for the PRBM model [27] of the Anderson localization. 
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